An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part II: Subcell finite volume shock capturing

The second paper of this series presents two robust entropy stable shock-capturing methods for discontinuous Galerkin spectral element (DGSEM) discretizations the compressible magneto-hydrodynamics (MHD) equations. Specifically, we use resistive GLM-MHD equations, which include a divergence cleaning mechanism that is based on generalized Lagrange multiplier (GLM). For continuous analysis to hold, and due divergence-free constraint magnetic field, system requires non-conservative terms, need special treatment. Hennemann et al. [DOI:10.1016/j.jcp.2020.109935] recently presented an strategy DGSEM Euler equations blends scheme with subcell first-order finite volume (FV) method. Our first contribution extension method systems such as In our approach, advective terms are discretized hybrid FV/DGSEM scheme, whereas visco-resistive only high-order We prove extended three-dimensional unstructured curvilinear meshes. derivation provides enhanced resolution by using reconstruction procedure carefully built ensure stability. provide numerical verification properties schemes meshes show their robustness accuracy common benchmark cases, Orszag-Tang vortex GEM reconnection challenge. Finally, simulate space physics application: interaction Jupiter's field plasma torus generated moon Io.